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The Ca 2+ -activated potassium channel KCa3.1 is emerging as a therapeutic target for a large variety of health disor- 
ders. One distinguishing feature of KCa3.1 is that the channel open probability at saturating Ca 2+ concentrations 
(Pomax) is low, typically 0.1-0.2 for KCa3.1 wild type. This observation argues for the binding of Ca 2+ to the 
calmodulin (CaM)-KCa3.1 complex, promoting the formation of a preopen closed-state configuration leading to 
channel opening. We have previously shown that the KCa3.1 active gate is most likely located at the level of the 
selectivity filter. As Ca 2+ -dependent gating of KCa3.1 originates from the binding of Ca to CaM in the C terminus, 
the hypothesis of a gate located at the level of the selectivity filter requires that the conformational change initiated 
in the C terminus be transmitted to the S5 and S6 transmembrane helices, with a resulting effect on the channel 
pore helix directly connected to the selectivity filter. A study was thus undertaken to determine to what extent the 
interactions between the channel pore helix with the S5 and S6 transmembrane segments contribute to KCa3.1 
gating. Molecular dynamics simulations first revealed that the largest contact area between the pore helix and the 
S5 plus S6 transmembrane helices involves residue F248 at the C-terminal end of the pore helix. Unitary current 
recordings next confirmed that modulating aromatic-aromatic interactions between F248 and W216 of the S5 
transmembrane helical segment and/ or perturbing the interactions between F248 and residues in S6 surrounding 
the glycine hinge G274 cause important changes in Pomax. This work thus provides the first evidence for a key 
contribution of the pore helix in setting Pomax by stabilizing the channel closed configuration through aromatic- 
aromatic interactions involving F248 of the pore helix. We propose that the interface pore helix/S5 constitutes a 
promising site for designing KCa3.1 potentiators. 



INTRODUCTION 



The voltage-insensitive calcium-activated potassium chan- 
nel of intermediate conductance, KCa3.1, has been doc- 
umented to play a prominent role in a large variety of 
physiological processes including immune reactions in- 
volving memory B and T cells (Wulff et al., 2004) , trans- 
epithelial ion transport in CP-secreting epithelial cells 
(Hayashi etal., 2012), control of vascular tone (Feletou 
and Vanhoutte, 2007), and proliferation/migration in 
various cell types (Cruse et al., 2006; Bradding and Wulff, 
2009; Gao et al., 2010). It follows that KCa3.1 is now 
recognized as a promising therapeutic target to treat 
life-threatening health disorders (Singh et al., 2001; 
Szkotak et al., 2004; Murthy et al., 2005; Roth et al., 2011). 
In line with this proposal are several reports demonstrat- 
ing that a pharmacological activation of KCa3.1 and KCa2.3 
channels constitutes a potential unique endothelium- 
specific antihypertensive therapy (Sheng et al., 2009; 
Damkjaer et al., 2012). Of interest is also the observa- 
tion thatKCa3.1 activation could improve cAMP-induced 
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CP secretion in tissues coming from cystic fibrosis 
patients with partial CFTR function, thus identifying 
KCa3.1 activation as a complementary strategy to correct 
the basic ion transport defect in cystic fibrosis (Roth 
etal., 2011). 

KCa3.1 channels are tetrameric membrane proteins 
with each subunit organized in six transmembrane 
segments, S1-S6, with a pore motif between segments 
5 (S5) and 6 (S6). The channel assembly and trafficking 
are regulated by the constitutively bound calmodulin 
(CaM) molecule, which also confers Ca 2+ sensitivity 
(Joiner et al., 2001). The crystal structure of the KCa3.1 
channel is not known. On the basis of homology model- 
ing data and substituted cysteine accessibility method 
(SCAM) experiments, we proposed that residues V275, 
T278, A279, V282, and A286 were lining the channel 
pore with residues C276, C277, L280, and L281, ori- 
ented opposite to the pore lumen (Simoes et al., 2002; 
Klein et al., 2007) . Additional experiments subsequendy 
demonstrated that the bundle-crossing region located 
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at the level of residues V282-A286 could not form a 
tight seal impermeable to K + ions for the closed KCa3.1 
channel. This conclusion came from experiments in which 
we showed that the accessibility of the thiol-modifying 
agent Ag + to cysteines engineered in the channel cav- 
ity was independent of the channel-conducting state 
(Garneau et al., 2009). In agreement with these results, 
Ba 2+ -based experiments on the Ca 2+ -activated KCa2.2 
channel demonstrated that the channel gate was most 
likely located deep in the channel central cavity (Bruening- 
Wright et al., 2007). In addition, data were presented 
indicating that negative gating of the KCa2.1 and 
KCa2.3 channels by NS8593 occurs through interac- 
tions with gating structures at a position close to the se- 
lectivity filter, deep in the pore inner vestibule (Jenkins 
et al., 2011). It was also suggested that the selectivity 
region of KCa2.x channels represents a binding site for 
potentiators such as GW542573X and CM-TPMF, with 
the residue equivalent to L215 in the KCa3.1 S5 trans- 
membrane helix playing a determinant role (Hougaard 
et al., 2009, 2012). Therefore, these observations tend 
to support a model where the selectivity filter region 
constitutes a key determinant to the action of KCa3.1 
regulators while demonstrating that residues distinct 
from the CaM-binding domain in the C terminus play a 
pivotal role in determining apparent Ca 2+ affinity and 
gating properties. However, channel gating is a dy- 
namic process that involves concerted conformational 
changes in many parts of the channel protein. As Ca 2+ - 
dependent gating of KCa3. 1 originates from the bind- 
ing of Ca 2+ to CaM in the C terminus, the hypothesis of 
a gate located at the level of the selectivity filter requires 
that the conformational change initiated in the C termi- 
nus be transmitted to the S6 transmembrane helix and 
to the channel pore helix directly connected to the se- 
lectivity filter. The exact molecular mechanism underly- 
ing KCa3.1 opening in response to Ca 2+ binding to the 
CaM-KCa3.1 complex in the C terminus remains to be 
elucidated. Structural information pertinent to channel 
gating was, however, obtained through the crystalliza- 
tion of CaM bound to the rat KCa2.2-CaM-binding do- 
main in the presence of Ca 2+ (Schumacher etal., 2001). 
On the basis of this structure, it was proposed that a large- 
scale conformational rearrangement is taking place in 
the presence of Ca 2+ , where the N-lobe of CaM binds 
to a C-terminal segment of an adjacent channel mono- 
mer resulting in a dimerization of contiguous subunits. 
This rearrangement would in turn lead to a rotation/ 
translation of the associated S6 transmembrane do- 
mains and to the opening of the ion-conducting pore 
(Schumacher et al., 2001, 2004; Wissmann et al., 2002; 
Maylie etal., 2004). 

As the CaM-binding domain of KCa3.1 is directly con- 
nected to the S6 transmembrane helix, activation of the 
channel gate at the level of the selectivity filter could 
depend upon the coupling between each of the channel 



pore helices and the associated S6 transmembrane seg- 
ment. SCAM experiments have already provided evi- 
dence for widespread conformational changes in the 
pore helix of the cyclic nucleotide-gated channel during 
gating (Liu and Siegelbaum, 2000). These observations 
were interpreted as evidence for a gate localized at the 
selectivity filter of the channel (Liu and Siegelbaum, 
2000). Mutations that perturb the pore helix of the 
Kir6.2 channel were similarly reported to affect the fast 
gating kinetics of the channel (Proks et al., 2001). The 
importance to gating of the interactions between the 
selectivity filter and the pore helix in channels was also 
confirmed in a study on KirBac3.1, which showed that 
changes at the level of the cytoplasmic domain were re- 
layed to the selectivity filter (Gupta et al., 2010). Simi- 
larly, changes in the conductivity of the KirBac selectivity 
filter were found to be correlated with remote confor- 
mational changes of the cytoplasmic domains, in agree- 
ment with the selectivity filter playing a role in gating by 
directiy responding to global conformational changes 
of the channel (Clarke et al., 2010) . Finally, work on the 
KcsA channel has led to the conclusion that C-type inac- 
tivation could arise from mechanical deformations prop- 
agating through a network of steric contacts between 
the channel inner helical bundle (equivalent to T278 
KCa3.1) and the C-terminal end of the pore helix (equiv- 
alent to L249-T250 KCa3.1) (Cordero-Morales et al., 
2007). Collectively, these observations suggest that resi- 
dues located deep in the KCa3.1 channel pore, close to 
the selectivity filter, could not only contribute to the 
CaM-KCa3. 1 complex to gating but also represent target 
sites for a pharmacological control of KCa3.1 activity. 

How the selectivity filter region in KCa3.1 is involved 
in gating remains unknown. However, one of the distin- 
guishing features of the KCa3.1 Ca 2+ dependence is that 
Pomax, the channel open probability at saturating Ca 2+ 
concentrations (>20 pM), remains low, typically 0.1-0.2 
for the wild-type channel, in contrast to a Pomax of 0.8 
reported for the KCa2.2 channel in the high activity 
mode (Hirschberg et al., 1998). The observation of a 
low Pomax at saturating Ca 2+ conditions argues for the 
binding of Ca 2+ to the CaM-KCa3.1 complex, promot- 
ing the formation of a closed-state configuration from 
which the channel transits to an open configuration. Such 
behavior has been documented in numerous ligand- 
gated channels where the agonist does not alter the 
open state but brings the channel to a preopen configu- 
ration. It is likely that the setting of Pomax depends on 
the energetics governing the interactions taking place 
at the selectivity filter level. 

In this study, we tested the hypothesis that interac- 
tions of the KCa3. 1 pore helix with the S5 and S6 trans- 
membrane segments contribute in setting Pomax. More 
specifically, we intended to identify which residues of 
the pore helix are involved in functionally coupling the 
pore helix to the S5 and S6 transmembrane segments 
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during gating. Our results demonstrate that important 
changes in Pomax can be obtained by modulating aro- 
matic-aromatic interactions between F248 of the pore 
helix and W216 of the S5 transmembrane helical seg- 
ment, and/ or by perturbing the interactions between 
F248 and the S6 transmembrane segment at the level of 
the G274 hinge residue. 

MATERIALS AND METHODS 

Computer-based homology modeling 

A first class of 3-D structures for the KCa3.1 pore region (F190- 
R287) was generated through homology modeling (Sali and 
Blundell, 1993) using as template the crystal structure of the Kvl.2 
channel (Protein Data Bank accession no. 2A79) . The choice of the 
Protein Data Bank 2A79 structure was essentially based on the analy- 
sis provided by the I-TASSER server, which identified Protein Data 
Bank 2A79 as the template with the best TM score of 0.869, know- 
ing that a TM score >0.5 corresponds approximately to two struc- 
tures of similar topology (Roy et al., 2010). The model includes 
part of the S4-S5 cytosolic linker plus the segment extending from 
the N terminus of the S5 transmembrane helix to the C terminus of 
the S6 transmembrane helix. A second class of models was ob- 
tained from the SAM-T08 server, which identified the Protein 
Data Bank accession number 3BEH (bacterial cyclic nucleotide- 
regulated channel) structure as the top-ranking template (Karplus, 
2009). Both templates were thus used to generate models of the 
KCa3.1 channel in the open (Protein Data Bank 2A79) and closed 
(Protein Data Bank 3BEH) configuration. 

Molecular dynamics (MD) simulations 

MD simulations were performed using the CHARMM-CGENFF 
force field with the channel pore region incorporated into a DPPC 
lipid bilayer in contact at top and bottom with an identical explicit 
water medium. Channel incorporation was accomplished accord- 
ing to the procedure implemented in the Charmm-Gui package 
(Jo et al., 2008). Overall, the system contains 55,732 atoms includ- 
ing 161 DPPC lipid molecules, 9,576 TIP3P model water molecules, 
plus 19 K + and 33 Cl~ ions to ensure electroneutrality at near physi- 
ological concentration. Cut-on and cutoff parameters needed to 
define nonbonded interactions were set to 10 and 12 A, respec- 
tively, and SHAKE constraints were used to fix lengths of bonds in- 
volving hydrogen atoms. The system was equilibrated according to 
the six-cycle scheme implemented in the Charmm-Gui output for a 
total equilibration period of 500 ps. Trajectories were generated for 
32-64 ns using a time step of 2 fs, and electrostatic and van der 
Waals interaction energies were computed from trajectories sam- 
pled at 0.2 ns. MD simulations were performed for a system at con- 
stant pressure (1 atm) and constant temperature (300°K). Solvent 
accessibility surface areas (SASAs) were averaged over a 64-ns tra- 
jectory using a standard probe of 1.6 A. 

Channel expression and molecular biology 

Wild-type and mutant KCa3.1 channels were expressed in Xenopus 
laevis oocytes by injection of their RNAs transcribed in vitro from 
pT7TS vector, which contains 5' and 3' untranslated regions of 
Xenopus [3-globin mRNA. Approximately 1-10 ng KCaS.l RNA was 
injected into each oocyte. Recordings were performed 1-7 d after 
injection. All the mutants were generated using the Q5 Site- 
Directed Mutagenesis kit (New England Biolabs, Inc.) and verified 
by sequencing at the Universite de Montreal genomics platform. 
Before patch clamping, the defolliculated oocyte was briefly incu- 
bated in a hyperosmotic solution containing (mM) 250 KC1, 
1 MgSC-4, 1 EGTA, 50 sucrose, and 10 HEPES, buffered at pH 7.4 



with KOH and the vitelline membrane peeled off using fine for- 
ceps. The oocyte was then transferred to a superfusion chamber 
for patch-clamp measurements. 

Solutions 

The bath and patch pipette solutions contained (mM) 200 K 2 S0 4 , 
1.8 MgCl 2 , 0.025 CaCl 2 , and 25 HEPES, buffered at pH 7.4 with 
KOH (referred to 200 mM K 2 S0 4 ). Sulfate salts were used to mini- 
mize the contamination from endogenous Ca 2+ -dependent chloride 
channels while enabling the chelation of contaminant divalent cat- 
ions such as Ba 2+ . Calcium-free solutions were prepared by adding 
1 mM EGTA to 200 mM K 2 S0 4 solutions without CaCl 2 . Bath solu- 
tion changes were performed as described previously using a rapid 
solution changer system (RSC-160; Biological) (Banderali et al., 
2004) . The solution exchange time was <10 ms. bis (2-mercaptoefhyl) 
sulfone (BMS; EMD Millipore) was directly dissolved in the 200-mM 
K 2 S0 4 solution at a concentration of 8 mM. 

Patch-clamp recordings 

Multiple- and single-channel inside-out recordings were per- 
formed using an amplifier (Axopatch 200A; Molecular Devices). 
Patch pipettes were pulled from borosilicate capillaries using a 
pipette puller (model PP-83; Narishige), which was used un- 
coated. The resistance of the patch electrodes ranged from 2 to 
5 Mfl. Data acquisition was performed using a Digidata 1320A 
acquisition system (Molecular Devices) at a sampling rate of 
1.0 kHz unless specified otherwise. Single-channel analysis was 
carried out using the QUB package (Qin et al., 1996, 1997). 
Dwell-time measurements were performed on data that were 
idealized according to the segmental-k means method based on 
a hidden Markov model-type analysis. 

Unnatural amino acids incorporation 

Incorporation of unnatural amino acids was performed accord- 
ing to the procedure described by Wang et al. (2006) . In brief, we 
used a pU6-pMpa plasmid (provided by P.G. Schultz, Scripps Re- 
search Institute, Lajolla, CA), which encodes two genes: Bacillus 
stearothermophilus Tyr-tRNA (with anti-codon mutated to CUA, 
complement to the TAG stop codon) and Escherichia coli Tyrosyl- 
tRNA-synthetase, which was mutated in Y37V/D182S/F183M to 
make the plasmid available to specifically recognize the unnatural 
amino acid jb-methoxy-L-phenylalanine (pMpa). The plasmid was 
coexpressed with the KCa3.1 channel mutated at position F248 
into an amber stop codon (TAG) in Xenopus oocytes, in which 5 mM 
pMpa (Sigma-Aldrich) had been injected. 

Statistical analysis 

Statistical significance was analyzed using unpaired Student's 
t test. P < 0.05 was considered statistically significant. Data are ex- 
pressed as mean ± SD. 

Online supplemental material 

Fig. SI shows a Western blot of full-length expression of a KCa3.1 
channel tagged with 6-His at position 132 in the external linker 
S1-S2, where residue F248 has been replaced by the unnatural 
amino acid pMpa. Fig. S2 presents a Western blot obtained for the 
F248W-W216F and W216F mutants, confirming channel expres- 
sion. Fig. S3 presents the open and closed dwell-time distributions 
computed for the wild-type, F248A, and F248W KCa3.1 channels. 
The online supplemental material is available at http://www.jgp 
.org/cgi/content/full/jgp.201311097/DCl. 

RESULTS 

Structural properties of the pore helix region 

The proposed structures for the 3BEH-based (bacterial 
cyclic nucleotide-regulated channel) and 2A79-based 
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(Kvl.2 channel) KCa3.1 channel pore region are pre- 
sented in Fig. 1 (A and B) . The selectivity filter is seen to 
be extending from T250 to T260 and is connected to 
the pore helix comprising amino acids G235 to L249. In 
both models, the residues V275, T278, V282, and A286 
(blue) are predicted to line the channel pore, in agree- 
ment with our previous SCAM studies (Klein et al., 2007) . 
An analysis of the distance between the Ca of corre- 
sponding residues along S6 on opposite subunits indi- 
cates that the narrowest part of the conduction pathway 
should be located at the level of V282, with a van der 
Waals pore diameter of ^4.5 A for the 3BEH-based 
model compared with >12 A for the model derived 
from the 2A79 template. It follows that the structure 
obtained from the 3BEH template is likely to be more 
representative of the KCa3.1 channel in the closed confor- 
mation relative to the 2A79-based model, which better 
accounts for the KCa3.1 open configuration. Both mod- 
els predict that residues 1244 and F248 of the pore helix 
should be projecting between the S5 and S6 transmem- 
brane helices, whereas L243 should be oriented toward 



the selectivity filter in proximity of V256 at the C-termi- 
nal end of the selectivity filter region. The proposed 
configuration for the pore helix was investigated in ex- 
periments where the spatial proximity of the L243 and 
V256 residues was tested through disulfide bond forma- 
tion between Cys engineered at L243 and V256, respec- 
tively. The inside-out patch-clamp recordings presented 
in Fig. 1 C confirmed in this regard that channel activ- 
ity of the L243C-V256C double mutant was sensitive to 
the addition of the small reducing agent BMS to the 
bathing solution. These observations confirmed that 
disulfide bond formation between L243C and V256C 
tends to stabilize the channel closed state configura- 
tion by constricting the selectivity filter. This effect was 
specific to the double mutation L243C-V256C system, as 
demonstrated in the multi-channel (L243C) and single- 
channel (V256C) recordings presented in Fig. 1 C 
where channel activity appeared insensitive to BMS ap- 
plication for single mutation channels. 

A computer analysis was next undertaken to deter- 
mine which residues of the KCa3.1 pore helix form the 
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Figure 1. Ribbon representation of the KCa3.1 pore 
region, including the S4-S5 cytosolic linker, the 
S5 and S6 transmembrane helices, the pore helix 
(gray), and the selectivity filter regions (stick repre- 
sentation). (A) Model structure generated by MOD- 
ELLERS. 11 using as template the crystal structure 
solved for a bacterial cyclic nucleotide-regulated 
channel MlotiKl (Protein Data Bank accession no. 
3BEH). Also shown is a space-filling representation 
of residues L243, F248, and V256, with the L243 side 
chain projecting opposite to the S5 and S6 trans- 
membrane segments toward V256 of the selectivity 
filter. In contrast, F248 is predicted to be oriented 
with its side chain making contact with both the S5 
and S6 helices. The residues documented to be fac- 
ing the channel pore are colored in blue (Klein 
et al., 2009) . (B) Model structure generated by MOD- 
ELLERS 11 using as template the crystal structure 
solved for the mammalian Kvl.2 channel (Protein 
Data Bank accession no. 2A79). As in A, residue 
F248 is seen as projecting between the transmem- 
brane helices S5 and S6. Only two subunits are pre- 
sented for clarity. (C) Inside-out patch recordings 
supporting the formation of a disulfide bond be- 
tween Cys engineered at 256 and 243. The reducing 
agent BMS caused channel activation when applied 
onto the L243C-V256C double mutation channel 
but failed to affect channel activity when applied 
to the V256C (single-channel) and L243C (multi- 
channel) single mutation channels. These observa- 
tions support the orientation of the pore helix relative 
to the selectivity filter proposed in A and B. Chan- 
nels were expressed in Xenopus oocytes, and current 
records were obtained at —60 mV in symmetrical 
200-mM K,S0 4 conditions at an internal Ca 2+ con- 
centration of 25 uM. The label "c" refers to the zero 
current level. Illustration by Discovery Studio Visual- 
izer (Accelrys). 
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interface with the S5 and S6 transmembrane segments. 
Our approach consisted of computing for each residue 
of the pore helix the difference in SASA (ASASA) with 
and without the S5 and S6 helices. It is expected that 
KCa3.1 residues not located at the interface formed by 
the pore helix with the S5 and S6 transmembrane seg- 
ments will have SASA values unaffected by the presence 
of the S5 or S6 helix (ASASA = 0) . In contrast, residues 
of the pore helix contributing to the interface with the 

55 and S6 helices should be characterized by a reduc- 
tion of their SASA values when calculations include 
either the S5 or S6 transmembrane segment. Fig. 2 
(A and B) summarizes the predicted contribution of 
each pore helix residue to the interface, with the S5 and 

56 transmembrane helices averaged over a 64-ns 
MD simulation. As seen, the largest contact areas with 
both the S5 and S6 helices involve residue F248 at the 
C-terminal end of the pore helix regardless of the tem- 
plate used. In contrast, residues T240, L241, and 1244 of 
the pore helix appeared to contribute more to the S5 
than S6 transmembrane segment interface, with ASASA 
values remaining overall smaller relative to F248. Notably, 



calculations for residues L249 and W242 revealed contact 
areas almost entirely limited to the S6 transmembrane 
helix, an effect coming from W242 and L249 contrib- 
uting to intersubunit contacts. As expected, ASASA 
values for L243 and 1246 predicted to be oriented fac- 
ing the selectivity filter were not affected by the presence 
of the S5 and S6 helices, indicating small if any inter- 
face contact surfaces. Collectively, this analysis showed 
that the largest contribution to the pore helix interface 
with both the S5 and S6 transmembrane segments is 
coming from residue F248. 

Because the CaM-KCa3.1 complex is direcdy con- 
nected to the S6 transmembrane helix, a similar proce- 
dure was applied to identify residues of the S6 segment 
contributing to the interface with the pore helix. The 
calculations presented in Fig. 2 (C and D) show that 
residues T278, C277, V275, G274, V272, G271, T270, 
and L268 in S6 are susceptible to interact with the chan- 
nel pore helix in both models. Interactions involving 
V275 and T268 appeared unique, however, as they rep- 
resent intersubunit contacts exclusively. Of interest, we 
noted that residue V275 can interface with L249 at the 
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Figure 2. Characterization 
of the interface formed by 
the pore helix with the S5 
and S6 transmembrane seg- 
ments based on the differ- 
ence in SASA (ASASA) for 
the KCa3.1 model derived 
from the 3BEH (A and C) or 
2A79 (B and D) templates. 
ASASA values were com- 
puted from MD 64-ns trajec- 
tories using the Charmm37 
force field. The contribution 
of individual pore helix resi- 
due to the interface formed 
by the pore helix and S5 
transmembrane segment or 
pore helix and S6 transmem- 
brane helix (A and B) was 
estimated as described in 
the Results. The results in A 
and B show that the average 
contact area with either the 

55 or S6 transmembrane seg- 
ment is maximal for residue 
F248 independently of the 
template used. In contrast, 
residues T240, L241, and 
1244 contributed more to the 
interface between the pore 
helix and the S5 than the 

56 segment, whereas residues 
W242 and L249 appeared to 

form an interface with S6 exclusively. Notably, residues L243 and 1246 failed to contribute to the interface formed by the pore helix with 
the S5 and S6 transmembrane segments, in accordance with these residues being oriented toward the channel selectivity filter. A similar 
procedure applied to the interface formed by the pore helix and the S6 transmembrane segment (C and D) revealed important contri- 
butions coming from T278, C277, V275, G274, V272, G271, T270, and L268, with the V275 and L268 contribution coming essentially 
from intersubunit interactions. It is concluded that the interface between the pore helix and the S5 plus S6 segments of each subunit 
essentially involves F248 at the C-terminal end of the pore helix. 
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Figure 3. Detailed represen- 
tations of the KCa3.1 channel 
pore helix region (gray) derived 
either from the 3BEH (bacterial 
nucleotide-activated potassium 
channel) or 2A79 (Kvl.2) tem- 
plates. A and B illustrate that 
the interface between F248 and 
the S5 transmembrane helix is 
template dependent, with F248 
predicted to be interacting with 
W216 and T212 of S5 for the 
3BEH-derived model but with 
L211 for the model obtained 
using the Kvl.2 structure (2A79) 
as template. C and D show in 
contrast a structural organiza- 
tion of the interface between 
F248 and the S6 transmem- 
brane segment in which F248 
is in both cases within proximity 
of the Gly hinge at 274 on S6 
added to potential contacts with 
C277, M273, and T270. Illustra- 
tions by Discovery Studio Visual- 
izer (Accelrys). 
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Figure 4. Electrostatic and 
van der Waals interaction 
energy between F248 and 
residues L268 to L280 of S6 
(A and B) or residues G210 
to T218 of S5 (C and D). En- 
ergy calculations were based 
on 64-ns MD trajectories ob- 
tained for the KCa3.1 model 
derived either from the 
3BEH (nucleotide-activated) 
or 2A79 (Kvl.2) templates. 
The bar graphs in A and B 
indicate that F248 of the pore 
helix strongly interacts with 
the Gly hinge at 274 for both 
models and to a lesser extent 
with T270, M273, C277, and 
T278. The interaction en- 
ergy pattern remained over- 
all template independent. In 
contrast, strong interaction 
energies were estimated be- 
tween F248 and W216 for the 
KCa3.1 model derived from 
the 3BEH template, whereas 
the maximum interaction 
energy for the KCa3.1 2A79- 
derived model involves resi- 
dues F248 and L211. This 
analysis thus indicates that the 
interaction energy pattern 
between F248 and residues 
in S5 varies according to the 
template used. 
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C-terminal end of the pore helix on an adjacent subunit. 
A close proximity between V275 and an adjacent pore 
helix would be in agreement with the observed strong 
channel inhibition by MTS reagents of the V275C mu- 
tant channel (Garneau et al., 2009), and with the con- 
tribution of V275 to the TRAM-34-blocking action 
(Wulff etal., 2000). 

Detailed representations of the interaction network 
between F248 and residues in S5 and S6 derived either 
from the 3BEH (bacterial nucleotide-activated potas- 
sium channel; A and C) or 2A79 (Kvl.2; B and D) tem- 
plates are presented in Fig. 3. In both models, the side 
chains of T270, G274, and C277 of the S6 transmem- 
brane helix are seen to be oriented toward the pore 
helix at the level of F248 (Fig. 3, C and D). However, the 
interaction pattern of the pore helix-S5 residues was 
found to be model dependent, with T212 and W216 
facing F248 in the 3BEH-based model (Fig. 3 A) , in con- 
trast to the structure derived from the Kvl.2 (2A79) 



template where F248 essentially makes contact with L21 1 
ofS5 (Fig. 3B). 

A detailed analysis of the nonbonded van der Waals 
plus electrostatic energy between F248 of the pore helix 
and residues extending from L268 to L280 along the S6 
transmembrane segment is presented in Fig. 4 (A and B) . 
The most important interactions were seen between 
F248 and the glycine hinge G274, arguing for F248 being 
within close proximity of the G274 backbone. Strong 
interactions were similarly estimated between F248 and 
T270, G271, M273, C277, and T278. These predictions 
appeared template independent and could be obtained 
for the 3BEH- and 2A79-based models. Interaction en- 
ergy between residue F248 and residues of the S5 trans- 
membrane helix is presented in Fig. 4 (C and D). The 
resulting interaction pattern appeared in this case 
model dependent, with F248 favorably interacting with 
W216, T212, and L215 in the 3BEH-derived structure, 
and with L211 for the Kvl.2-based open configuration. 
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Figure 5. Alanine scan of the KCa3.1 pore 
helix region. (A) Inside-out current record- 
ings obtained in symmetrical 200-mM K 2 S0 4 
conditions at a pipette potential of 100 mV 
(T240A, I244A, and F248A) or 150 mV 
(T250A) from KCa3.1 mutants expressed in 
Xenopus oocytes. The label "c" refers to zero 
current level. The internal Ca 2+ concentra- 
tion was 25 pM throughout. The mutation 
F248A caused a highly significant increase of 
the channel open probability at saturating in- 
ternal Ca 2+ concentration (Pomax) relative to 
wild type (P < 0.0005) compared with T240A 
(P < 0.01). The substitution I244A caused a 
nonsignificant variation in Pomax relative to 
wild type. Because of its strategic position at 
the N-terminal end of the selectivity filter, the 
mutation of T250 to Ala resulted in an impor- 
tant decrease in channel conductance (9 pS 
compared with 30 pS), with Pomax being not 
significantly different from wild type. These 
results are summarized in B with Pomax val- 
ues of 0.42 ± 0.08 (n = 4), 0.10 ± 0.04, 0.75 ± 
0.01 (n = 3), and 0.16 ± 0.01 (n = 2) for the 
T240A, I244A, F248A, and T250A channels, 
respectively. Our results point toward the 
interactions involving F248 with the S5 and 
S6 helices as being determinant in setting 
Pomax. Also shown are all-point histograms 
computed from the entire recording. 
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Which residues of the pore helix affect the channel 
maximum open probability? 

The previous analysis has led to the identification of 
residues predicted to line the interface formed by the 
pore helix with the S5 and S6 transmembrane segments. 
An alanine scan was thus performed to determine how 
important these residues are in setting Pomax. Typical 
single-channel recordings are presented in Fig. 5 A. We 
note that the mutation F248A resulted in a Pomax of 
0.75 ±0.10 (w=3) compared with 0.22 ± 0.07 (w=8) for 
the wild-type channel (Fig. 5 B, WT). The mutation 
T240A led to a less drastic increase in Pomax with a 
value of 0.42 ± 0.08 (n= 4) compared with I244A, which 
resulted in a Pomax decrease to 0.10 ± 0.04 (n = 4). The 
mutation T250A had a minor impact on the channel 
Pomax but resulted in an important decrease in the 
channel unitary conductance, from 30 pS for KCa3.1 
wild type to 9 pS for the T250A mutant (not depicted) . 



Collectively, these results show that the mutation F248A 
is unique, as it caused a 3.2-fold increase in Pomax. 

Aromatic residues control Pomax 

To determine if hydrophobicity and/ or volume of the 
residue at position 248 control Pomax, experiments 
were performed in which F248 was substituted by resi- 
dues differing in size and/ or hydrophobicity. Typical 
unitary current traces are presented in Fig. 6. These 
recordings demonstrate that among all the mutants 
tested, F248T only yielded Pomax values (0.19 ± 0.05; 
n= 3) not superior to wild type (F248, 0.22 ± 0.07; n= 8) . 
Clearly, there was no direct correlation between Pomax 
and the volume of the substituting residue, as replacing 
F248 by a smaller (F248A) or larger amino acid (F248W) 
resulted in both cases in an increased Pomax relative to 
wild type. In addition, we noted that although Val (V) 
and Thr (T) are of roughly the same shape and volume, 




Figure 6. Examples of inside-out current recordings (A) and calculated Pomax (B) for KCa3.1 channels after the substitution of F248 
by residues differing in volume size and/or hydrophobicity. Recordings were performed in symmetrical 200-mM K 2 S0 4 conditions in 
saturating internal Ca 2+ (25 pM) at a pipette potential of 100 mV. Low Pomax values compared with wild type were obtained with the 
F248T mutant only with a mean value of 0.19 ± 0.05 (n= 3). The substitution of F248 by either a smaller (F248A) or larger (F248W) resi- 
due led systematically to a Pomax increase. Also shown is an example of inside-out current trace obtained in conditions where F248 was 
replaced by the unnatural amino acid pMpa, a tyrosine analogue obtained by substituting the hydroxyl moiety of Tyr by an 0-CH 3 group. 
This substitution led to a drastic increase in Pomax with a mean value of 0.94 ± 0.03 (n = 4), compared with 0.22 ± 0.07 (n = 8) for wild 
type. (B) Bar graph summarizing the effect of substituting F248 by nonpolar residues differing in volume size and/or hydrophobicity. 
The resulting Pomax ranking reads: T~F<H<Y<S<A<V<L~W< pMpa, confirming the absence of correlation between Pomax 
and the volume size and/or hydrophobicity (A<S<T<H<V<Y<L<F< W) (Wimley et al., 1996) of the substituting residue. Also 
shown are all-point histograms computed from the entire recording. Bars, 2 pA and 0.5 s. 
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there was a significant difference in channel activity, 
with Pomaxof0.77± 0.05 (n = 3) and 0.19 ±0.05 (w=3) 
for the F248V and F248T mutants, respectively. We 
concluded that the substitution of a methyl (V) by a hy- 
droxyl group (T) modified the interaction energy pat- 
tern involving residue 248 as to cause a decrease in 
Pomax value. The results of this analysis are summa- 
rized in Fig. 6 B, which shows that the overall Pomax 
ranking reads T = F<H<Y<S<A<V<L=W. 

A scatter plot of Pomax as a function of residue 248 
volume size is presented in Fig. 7 A. This analysis reveals 
two distinct sets of data. Whereas there was a modest 
but significant increase in Pomax as a function of vol- 
ume size for Ala, Val, and Leu, drastic changes were ob- 
served with the aromatic (Phe, Tyr, His, and Trp) residues. 
Clearly, increasing the size of the aromatic side chain 
engineered at position 248 resulted in higher Pomax 
values. Importantly, aromatic residues such as Phe and 
Tyr showed significant lower Pomax values compared 
with nonaromatic amino acids of similar sizes such as Val 
or Leu, suggesting specific effects related to the pres- 
ence of aromatic rings in these cases. To confirm that 
the size of aromatic residues at 248 is a key determinant 
to Pomax, experiments were undertaken in which we used 
the unnatural amino acid pMpa, a tyrosine analogue 
obtained by substituting the hydroxyl moiety of Tyr by 
a 0-CH 3 group. The resulting aromatic amino acid has 
an estimated volume of 227 A 3 compared with 192 A J 
for Tyr and 203 A' ! for Trp. Single-channel recordings 
confirmed that the substitution of F248pMpa caused an 
important increase in Pomax, with a mean value of 0.94 + 
0.03 (n = 4), in agreement with the proposal of a strong 
correlation between Pomax and the volume of the aro- 
matic residue engineered at 248 (see Fig. 7 A). 



Figure 7. (A) Effect on Pomax of substituting F248 by either 
nonaromatic or aromatic residues. Scatter plot of Pomax as a 
function of the residue volume. This analysis reveals two distinct 
sets of data. Whereas there was a modest but significant increase 
in Pomax as a function of volume size for Ala (A) , Val (V) , and 
Leu (L), drastic changes were observed with the aromatic (Phe 
[F], Tyr [Y], His [H], Trp [W], and pMpa) residues. Importantly, 
aromatic residues such as Phe (F) and Tyr (Y) showed significant 
lower Pomax values compared with nonaromatic amino acids 
of similar sizes, such as Val (V) or Leu (L), suggesting specific 
effects related to the presence of an aromatic residue at 248. 
(B and C) Ln(Keq_mutant/Keq_wild_type), where Keq = Pomax/ 
(1 — Pomax) plotted as a function of the energy perturbation 
AAEoc coming from mutating residue 248 for the F248T, F248V, 
F248A, and F248W mutants. AAEoc was calculated as AEoc_ 
mutant — AEoc_wild_type, where AEoc is the difference in non- 
bonded energy for residue 248 between the open and closed 
state. Nonbonded energies were computed from >32-64-ns MD 
trajectories by averaging the van der Waals plus electrostatic in- 
teractions between the residue at 248 and the surrounding atoms 
within a 10-A radius (excluding lipids and water) (B), or between 
the residue at 248 and the S5 and S6 helices (C) . The resulting 
analysis is in agreement with the Pomax ranking F248T = F248 < 
F248A b F248V < F248W observed experimentally, while estab- 
lishing an energy equivalence (AAEoc « 0) between a Thr (T) 



S5 transmembrane helix control of Pomax via F248 

The observation of a lower Pomax value for aromatic 
compared with nonaromatic residues engineered at 
248 could be indicative of specific interactions involving 
aromatic-aromatic interactions potentially contribut- 
ing to the energy barrier controlling channel opening. 
Aromatic amino acids are known to play a prominent 
role in determining protein structure through tt-tt or 
hydrogen-Tr interactions. For polar hydrogen-ir inter- 
actions, the donors are hydrogen atoms connected to 



and a Phe (F) at 248. The latter result argues for tt-tt interac- 
tions involving F248 and the S5 transmembrane segment in the 
closed state as being determinant in setting Pomax for the KCa3. 1 
wild-type channel. These simulations also predict an increase in 
Pomax relative to wild type with the F248W mutant despite the 
possibility of tt-tt interactions with S5, an effect attributable to 
the mutation F248W stabilizing to a greater extent the channel 
open than closed configurations. In contrast, the Pomax increase 
seen with the F248A and F248V mutants appeared to be related to 
a stronger destabilization for the channel closed than open state. 



Garneau et al. 297 



electronegative atoms (for instance, R-OH, R-NH>) , and 
the acceptors may be represented by various aromatic 
molecules and conjugate it groups (Du et al., 2013). 
More importandy, we note that the 3BEH-based model 
of KCa3.1 illustrated in Fig. 3 A shows possible face- 
edge interactions between F248 and W216 of the S5 
transmembrane helix. As Tyr, Tip, His, and Phe constitute 
polar hydrogen acceptors, which can potentially interact 
with a polar hydrogen donor such as Trp (NH) (Du et al., 
2013), the effect on Pomax seen with the substitutions 
F248Y/H/W could be indicative of polar hydrogen-n 
interactions and/ or tt-tt interactions between the resi- 
due 248 in the pore helix and W216 in S5. Such interac- 
tions would be absent with the F248A/V/L mutations, 
leading to higher Pomax values as seen experimentally. 
As Trp can also behave as a hydrogen acceptor in polar 
hydrogen-ir interactions, a strong interaction between 



F248T (hydrogen donor) and W216 could explain the 
unexpected observation of F248T leading to a smaller 
Pomax compared with F248A despite a larger side-chain 
volume. However, this interaction would be size depen- 
dent, as the Pomax for the F248S mutant is similar to 
the Pomax of F248A. It follows that substituting W216 
(hydrogen donor) by nonaromatic residues such as Ala 
or Leu should impair aromatic-aromatic interactions 
involving F248 and reproduce the effects observed when 
F248 was substituted by a nonaromatic amino acid. 
Examples of single-channel recordings obtained with 
the W216A and W216L mutants are illustrated in 
Fig. 8 A. Clearly, replacing W216 by a nonaromatic resi- 
due caused a drastic increase in Pomax, with values of 
0.68 + 0.04 (n = 3) and 0.95 ±0.01 ( n = 3) for W216A 
and W216L, respectively (see Fig. 8 C). The effect ob- 
served with W216A cannot be attributed to a decrease 
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Figure 8. Evidence indicating potential aromatic-aromatic interactions involving residue 248 of the pore helix and W216 of S5. 
(A) Inside-out current recordings demonstrating that the absence of aromatic-aromatic interactions resulting from the substitution 
of W216 by the nonaromatic Ala or Leu residue caused a strong increase in channel activity, with Pomax values of 0.68 ± 0.04 (re = 3) 
and 0.95 ± 0.01 (n = 3), respectively, compared with 0.22 ± 0.07 (re = 8) for the wild-type channel (W216). However, such behavior was 
not observed with the L215A mutant nor the T212A (not depicted). Substituting W216 by the smaller Phe aromatic residue (W216F) 
resulted in a single-channel fluctuation pattern characterized by a period of high activity (Pomax = 0.70 ± 0.2; re = 3), ending within 
<3 min in 50% of the recordings ( re = 10) by a transition to stable closed state. Such behavior is compatible with strong aromatic-aromatic 
F248-W216F interactions maintaining the channel in the closed configuration. These results are summarized in the bar graph presented 
in C. These observations agree with a model whereby aromatic-aromatic interactions between F248 of the pore helix and W216 of S5 
tend to stabilize the channel closed configuration, as suggested by the 3BEH-based model illustrated in B. Recordings were performed 
in symmetrical 200-mM K 2 S0 4 conditions in saturating internal Ca 2+ (25 jiM) at a pipette potential of 100 mV. All-point histograms were 
computed from the entire recording. Illustration by Discovery Studio Visualizer (Accelrys). 
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of the side-chain volume, as the substitution W216L led 
to a higher Pomax value relative to W216A. Collectively, 
these results support a model in which polar hydrogen-ir 
and/or tt-tt interactions involving residue 248 of the 
pore helix and W216 on S5 (Fig. 8 B) contribute to sta- 
bilize the gate closed configuration, while confirming 
the prominent role of the selectivity filter region to 
the KCa3.1 gating process. In line with this proposal is the 
observation that substituting W216 by the smaller Phe 
aromatic residue (W216F) resulted either in the ab- 
sence of detectable channel activity (50% of the 25 oo- 
cytes tested) despite Western blot evidence for channel 
expression (see Fig. S2), or in single-channel records 
where the channel gate switched within <3 min in 50% of 
the recordings considered, from a high Pomax state (0.7 ± 
0.2; n=S) to a permanent closed configuration (Fig. 8 A) . 
Such behavior was not seen when W216 was mutated to 
a nonaromatic residue such as Ala or Leu (Fig. 8 A) . These 
observations are compatible with aromatic-aromatic in- 
teractions involving F248 and W216F significantly in- 
creasing the likelihood of forming a stable closed state. 



Notably, experiments in which the Phe at 248 and Trp at 
216 were interchanged (F248W-W216F) failed to yield 
detectable channel activity despite biochemical evi- 
dence of channel expression (see Fig. S2). 

Other residues in S5 could be interacting with F248 

The structural model and the interaction energy dia- 
gram illustrated in Figs. 8 B and 4 C, respectively, point 
toward T212 and L215 in S5 as potential residues sus- 
ceptible to be within the interaction range of F248. The 
contribution of these two residues in setting Pomax was 
investigated in single-channel experiments using the 
T212A and L215A mutant channels. In contrast to the 
mutation W216A, the substitution T212Aor L215Afailed 
to significantly affect channel activity, with Pomax val- 
ues of 0.28 ± 0.15 (w=3) and0.19±0.04 (n= 3), respec- 
tively (Fig. 8 C). These observations do not support a 
model whereby OH-tt interactions involving T212 and 
F248 would participate in the control of Pomax. In con- 
trast to the model derived from the 3BEH template, the 
model based on the Kvl.2 (2A79) structure predicts 
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Figure 9. Effects of mutating residue L211 in S5 on KCa3.1 Pomax. (A) Inside-out current recordings demonstrating that increasing 
the volume size of the substituting residue at position 211 in S5 has no significant impact on Pomax relative to wild type. A significant 
increase was, however, observed with the L211F, suggesting possible aromatic-aromatic interactions between L211F and F248. The 
results of this analysis are presented in the scatter plot illustrated in B, which shows a Pearson correlation coefficient of 0.75 ± 0.15 
between the volume size of the residue engineered at 211 and Pomax. (C) A better correlation was obtained between Pomax and the 
solvation energy of the residue at 211 with a Pearson correlation coefficient of 0.84 ± 0.07. This analysis confirms that L211 can come in 
close proximity of F248 during gating, as revealed by the Pomax increase seen with the L211F mutant. Recordings were performed in 
symmetrical 200-mM K 2 S0 4 conditions in saturating internal Ca 2+ (25 pM) at a pipette potential of 100 mV. All-point histograms were 
computed from the entire recording. 
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strong interaction energy between F248 and L211, and 
to a lesser extent with L215. The importance of these 
interactions in setting Pomax was investigated in inside- 
out patch-clamp experiments where L21 1 was substituted 
by residues of variable volume size and hydrophobicity. 
The results of these experiments are summarized in Fig. 9. 
On the basis of the model derived from the Kvl.2 tem- 
plate, one would expect a potential stabilization of the 
channel open state by favoring tt-tt or OH-tt interac- 
tions between F248 and Phe (F), Trp (W), or Thr (T) 
residues engineered at 211. Our results indicate a small 
increase in Pomax with the L211F mutant relative to 
wild type, whereas the mutations L211T and L211W ap- 
peared ineffective in modifying Pomax. At best, we found 
a correlation between the volume size of the residue en- 
gineered at 211 and Pomax, with a Pearson correlation 
coefficient of 0.75 ± 0.15 (Fig. 9 B). A stronger correla- 
tion was obtained between Pomax and the solvation en- 
ergy of the residue engineered at 211 with a Pearson 
coefficient of 0.84 ± 0.07 (Wimley et al, 1996) (Fig. 9 C) . 
Collectively, these observations suggest that the interac- 
tions between F248 and the S5 transmembrane helix 
are essentially governed by W216 for the channel in the 
closed configuration (3BEH-based model) , with no sig- 
nificant contribution from L211, T212, and L215. 

F248 interactions with the S6 transmembrane helix 
The results presented in Fig. 4 indicate that F248 can 
interact preferentially with the Gly hinge at 274 and 
with residues T270, M273, C277, and T278 of the S6 
transmembrane helix. MD simulations further suggest 
that the strong interactions with G274 arise in part from 
van der Waals contacts indicating a very close proximity 
between the G274 backbone and the F248 side chain. 
As the CaM-KCa3.1 complex is directly connected to 
the S6 transmembrane segment, a double mutant cycle 
analysis was performed to determine how these interac- 
tions could contribute to Pomax. In this analysis, the cou- 
pling energy AAG was computed as (AG mutati(>n i_ mutatloll2 + 

AGwT-Wt) ~~ (AGmutationl-WT + AG mu tation2-Wr) ! with AG 

given by kT In (Pomax) , with k the Boltzmann constant 
and T the temperature (Horovitz, 1996; Wall-Lacelle 
et al., 2011). Examples of cycle analyses are shown in 
Fig. 10. We note first the absence of coupling with ei- 
ther F248 or F248W for the C277A, M273A, and T270A 
mutants with AAG of 0.02 ± 0.28 kcal/mol, -0.27 ± 0.30 
kcal/mol, and -0.01 ± 0.10 kcal/mol for the F248W- 
C277A, F248W-M273A, and F248W-T270A double 
mutant cycle systems, respectively (not depicted) . In ad- 
dition, the observation that the mutations T270A and 
T270S were ineffective in modifying Pomax when added 
to the F248 or F248W template (not depicted) was seen 
as evidence for T270 not being involved in OH-iT-type 
interactions with F248 or F248W. Strong AAG of - 1 .04 ± 
0.26 kcal/mol and —1.15 ± 0.56 kcal/mol was obtained, 
however, with the F248W-M273W and F248W-C277W 



double mutants (Fig. 10). No current could be detected 
with the T270W mutant (not depicted) . Of interest, the 
mutation F248W failed to cause an increase in Pomax 
when added to the C277W template, an effect not seen 
with C277A (not depicted). These observations sug- 
gest potential aromatic-aromatic interactions between 
F248W and C277W, as seen with W216 (see Fig. 1 1) . We 
also note that the mutation M273W (Fig. 10 A) led to a 
Pomax of 0.95 ± 0.005 (n = 3) with F248 and to 0.63 ± 
0.05 ( n = 4) when F248 was substituted by a Trp (F248W- 
M273W). With a Pomax of 0.86 ± 0.07 (n = 5) for 
F248W-M273, 0.80 ± 0.10 (n = 3) for F248W-M273A 
(not depicted), and 0.63 ± 0.05 (n = 4) for F248W- 
M273W (Fig. 10 B), there is no apparent correlation 
between the volume size of residue 273 and Pomax in 
conditions where F248 has been mutated into a Trp. 
Such behavior was not found when mutations of M273 
were performed on the wild-type KCa3.1 template (F248) 
with a Pomax of 0.22 ± 0.02 (n = 3) and 0.95 + 0.005 
(n = 3) for the F248-M273A (not depicted) and F248- 
M273W (Fig. 10 B) channels, respectively. Clearly, the 
presence at 248 of a Phe with a smaller side-chain vol- 
ume compared with Trp caused the channel Pomax to 
become more sensitive to the volume of the residue sub- 
stituting for M273. These results confirm the close 
proximity between F248 and M273, as predicted from 
the homology-based representations of the KCa3.1 open 
and closed configurations (see Fig. 11). 

DISCUSSION 

The results presented here demonstrate that aromatic- 
aromatic interactions between F248 of the pore helix 
and W216 of the S5 transmembrane segment play a cru- 
cial role in setting the channel maximum open proba- 
bility, while supporting the proposal of an active gate 
located at the level of the selectivity filter. It is also shown 
that F248 is within the interaction range of C277 and 
M273 adjacent to the Gly hinge of the S6 transmem- 
brane segment and could thus contribute to the func- 
tional link between the CaM-KCa3.1 complex and the 
pore helix. These results suggest that the channel closed 
configuration is stabilized by aromatic-aromatic inter- 
actions between F248 at the C-terminal end of the pore 
helix and the S5 transmembrane segment, and support 
a model where the interface between the pore helix 
and the S5 transmembrane segment constitutes a po- 
tential target site for the design of KCa3.1 potentiators. 

Modeling of the KCa3.1 pore region 

Our homology modeling analysis has led to the iden- 
tification of two templates susceptible to account for 
the KCa3.1 pore region: the mammalian Shaker Kvl.2 
channel (2A79) (Long et al., 2005a) and the bacterial 
cyclic nucleotide-regulated MlotiKl channel (3BEH) 
(Clayton et al., 2008). The resulting structures likely rep- 
resent two distinct configurations of the KCa3.1 channel 
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with the 2A79- and 3BEH-based structures referring to 
the channel open and closed states, respectively. SASA 
calculations have also revealed that the largest contact 
areas with both the S5 and S6 helices involve residue 
F248 at the C-terminal end of the pore helix regard- 
less of the template used. The results of these analyses 
strongly suggest that F248 participates in conveying the 
structural changes induced by the binding of Ca 2+ to the 
CaM-KCa3.1 complex to the pore helix. This proposal 
is in line with data coming from mass spectroscopy mea- 
surements where it was demonstrated that residue M94 
of KirBac3.1 (equivalent to F248 in KCa3.1) displays 
among the largest changes in solvent accessibility when 
changing from a closed to an open configuration (Gupta 
et al., 2010). These observations were taken as evidence 
for a channel gate located close to the selectivity filter 
(Gupta et al., 2010), in accordance with the present re- 
sults demonstrating a key contribution of F248 in setting 
the channel Pomax. Changes in solvent accessibility 



during KCa3.1 gating are also expected on the basis of 
our SASA calculations, with the surface contribution of 
1244 and F248 to the interface formed by the pore helix 
with the S5 and S6 transmembrane segments varying 
from 114 and 173 A 2 for the closed 3BEH-based config- 
uration to 102 and 162 A 2 for the open-like configuration 
(unpublished data). In contrast, the surface contribu- 
tion of L249 did not change appreciably between the 
two channel model structures, with respective values of 
99.4 and 99.2 A 2 (unpublished data). Thus, these calcu- 
lations support an increase in solvent accessibility at 1244 
and F248 when changing from a closed to an open state, 
an effect not seen at L249, which is more exposed to the 
channel central cavity. 

Aromatic interactions between the pore helix and 
the S5 transmembrane segment 

Whereas the 3BEH- and 2A79-based models led to simi- 
lar interaction energy profiles between F248 and S6 
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Figure 10. Mutant cycle analyses demonstrating functional coupling between F248W of the pore helix and C277W or M273W of the 
S6 transmembrane segment. Energies are expressed in kilocalorie/mole. (A) Example of inside-out current recordings illustrating the 
effect of the M273W mutation on the change in Pomax resulting from the substitution at 248 of a Phe by a Trp. These observations are 
summarized in the cycle diagram illustrated in B, which shows that the mutation M273W resulted in a drastic Pomax increase (0.95 ± 
0.005; n = 3) , an effect that was partly impaired with a Trp at 248 (Pomax of 0.63 ± 0.05; ra=4). (C) Examples of inside-out current record- 
ings illustrating the coupling between Trp at 248 and 277, respectively. The cycle diagram in D indicates in this regard that the muta- 
tion C277W totally prevented the Pomax increase normally resulting from the substitution F248W, while being ineffective in modifying 
Pomax with F248. These observations would be compatible with Trp-Trp interactions between F248W at the pore helix and C277W on 
S6 stabilizing the closed state (Fig. 11). On the basis of this analysis, it is concluded that both C277 and M273 are in proximity of F248. 
Recordings were performed in symmetrical 200-mM K 2 S0 4 conditions in saturating internal Ca 2+ (25 uM) at a pipette potential of 
100 mV. All-point histograms were computed from the entire recording. 
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residues, important differences were seen in the inter- 
action energy profiles between F248 and residues of the 
S5 transmembrane helix. Of interest is the observation 
that F248 is mostly interacting with W216 and T212 of 
S5 for the 3BEH-based configuration, in contrast to the 
2A79-based structure where F248 is predicted to inter- 
act for the most part with L21 1 and L215. More impor- 
tantly, our results support a model whereby interactions 
between aromatic residues at 248 and W216 control 
Pomax. Substituting either F248 or W216 by nonaro- 
matic amino acids systematically led to an increase in 
Pomax, arguing for aromatic-aromatic interactions sta- 
bilizing the channel closed configuration. Aromatic 
rings have been documented to interact preferably in a 
stacked offset mode or T-shape mode in which the H 
atoms on one ring points toward the ir-electron cloud of 
the second ring (Samanta et al, 1999) . Our results would 
be in accordance with tt-tt and/ or polar hydrogen- 
Tr-type interactions accounting for Phe, His, and Tyr at 
position 248 in the pore helix showing significant lower 
Pomax values compared with nonaromatic amino acids of 
similar sizes, such as Val or Leu. These aromatic-aromatic 
interactions are likely to stabilize the selectivity filter 
closed configuration leading to low Pomax values. A de- 
crease in Pomax relative to wild type was also observed 
with the F248T mutant. On the basis of the side-chain 
volume alone, the F248T mutation should have led to a 
Pomax comparable to F248V, which is clearly not the 
case (see Fig. 6) . Additional forces might contribute in 
this case to the interactions with Trp; in particular, a 
hydroxyl group close to the tt face of W216 could be 
responsible for polar hydrogen-TT interactions. It is also 
apparent from the data presented in Figs. 7 A and 8 that 
increasing the volume size of the residues at 248 re- 
sulted in higher Pomax values, an effect that was par- 
ticularly strong when coupled to interactions involving 
aromatic rings. 

To get further insights on how nonbonded interac- 
tions involving F248 control Pomax, 32-64-ns MD simu- 
lations were performed for the wild-type (F248) , F248T, 
F248A, F248V, and F248W mutant channels in the closed 
(3BEH template) and open (2A79 template) configura- 
tions. The resulting trajectories were used to compute 
the energy perturbation AAEoc coming from mutating 
residue 248, where AAEoc = AEoc_mutant-AEoc_wild_ 
type, with AEoc as the difference in nonbonded energy 
for residue 248 between the channel open and closed 
states. AEoc values were obtained by averaging the non- 
bonded interaction energy (van der Waals and electro- 
static) between residue 248 and the surrounding atoms 
within a 10-A radius (excluding lipid and water) (Fig. 7 B) 
or between residue 248 and the S5 plus S6 transmem- 
brane helices exclusively (Fig. 7 C) . It follows that the 
ratio Keq_mutant/Keq_wild type, where Keq = Pomax/ 
(1 — Pomax) (Fig. 7 A) should vary according to 
exp(— AAEoc/kT), as long as mutation at 248 does not 



affect to a significant degree the nonbonded interac- 
tions between residues other than 248. Furthermore, 
AAEoc values refer to energy changes for a single sub- 
unit and do not include a possible contribution coming 
from a combination of subunits. Therefore, negative 
AAEoc values reflect an increase in Pomax relative to 
wild type as the ratio Keq_mutant/Keq_wild type be- 
comes >1 under these conditions. The results of these 
calculations are illustrated in Fig. 7 (B and C), where 
the logarithm of the ratio Keq_mutant/Keq_wild type is 
plotted as a function of AAEoc. Despite limitations com- 
ing from the use of homology-derived structures de- 
scribing the selectivity filter region only, our trajectory 
analysis demonstrates a clear correlation between AAE oc 
and the logarithm of the ratio Keq_mutant/Keq_wild_ 
type, in agreement with our experimental observation 
of a Pomax ranking corresponding to F248 ~ F248T « 
F248A « F248V < F248W. This correlation holds for 
AAE oc computed for residue 248 interacting with all the 
atoms within a 10-A radius and for 248 interacting with 
the S5 and S6 helices exclusively. Furthermore, these 
calculations indicated that in contrast to Ala (A) and 
Val (V), the mutation F248T caused an equal stabiliza- 
tion of the channel closed (E dosedF248T — E dosedF248 = 
— 7.91 kcal/ mol, with E dosed mutant = van der Waals + elec- 
trostatic interaction energy between residue 248 and 
atoms in a 10-A radius for the closed configuration) and 
open state (E openF248T - E openF248 = -7.96 kcal/ mol) , in 
accordance with AAE oc ~ 0 (Fig. 7 B) in this case and a 
Pomax (F248T) ~ Pomax (wild type) as observed exper- 
imentally (Fig. 7 A) . The absence of a significant energy 
difference AAE oc with F248T and not with the F248A 
and F248V mutants supports of a model whereby aro- 
matic-aromatic interactions involving F248 are deter- 
minant in setting Pomax for the KCa3.1 wild-type 
channel. We note, however, that the mutation F248S 
has led to a Pomax value comparable to F248A despite 
the possibility of OH-tt interaction as for F248T. MD 
simulations confirmed in this regard that the F248S mu- 
tant behaves as F248T and should be characterized by a 
low Pomax, as observed with F248T. However, the mu- 
tation F248S has appeared unique among the mutant 
channels tested, including F248T, as it not only affected 
the channel Pomax but also caused a major decrease 
(threefold) of the channel unitary conductance (10 pS), 
arguing for a direct effect on the selectivity filter. Simi- 
larly, the F248C mutant where the O-H group of Ser is 
now replaced by S-H failed to yield detectable channel 
activity (not depicted), suggesting again a robust im- 
pact on the selectivity filter region. We speculate that 
these additional interactions and how they can be re- 
lated to the channel Pomax cannot be properly taken 
into account by our current models, thus limiting the 
predictive value of MD simulations for the F248S mu- 
tant. The results in Fig. 7 (B and C) equally show that 
the energy difference AAE oc for the F248W channel is 
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within the energy range computed for the nonaromatic 
Ala (A) or Val (V) residues despite the possibility in this 
case of it-it interactions as previously discussed. An 
analysis of the associated MD trajectories revealed that 
the mutation F248W better stabilized the channel open 
(E open F248W ~~ E op en_F248 = — 7.5 kcal/mol) than the closed 
state (E closed F2 48w - E dosed _ F248 = -4.0 kcal/mol), result- 
ing in AAE oc = —3.4 kcal/mol for an energy balance 
favoring an increase in Pomax relative to wild type as 
observed experimentally. This effect could be linked to 
the mutation F248W leading to an increase in stabiliz- 
ing van der Waals interactions with the S5 and S6 heli- 
ces in the open state, a mechanism that would account 
for the results presented in Fig. 7 A, where Pomax is 
seen to increase as a function of the residue volume for 
aromatic amino acids. A stabilization of the open state 
would also be compatible with the conclusions coming 
from a dwell-time analysis that confirmed that the muta- 
tion F248W results in a drastic right displacement of the 
open-time distribution curve from a mean value of 7 ms 
for wild type to 285 ms for F248W (Fig. S3). Changes 
in the closed-time distributions were equally observed 
characterized by the appearance of a dominant short- 
time distribution centered at 3 ms in support of a main 
destabilization of the channel closed configuration 
(Fig. S3). The Pomax augmentation observed with the 
F248A mutant seems, however, to arise from a different 
mechanism involving a greater destabilization of the 
channel closed state (E dosed _ f248 a — E dosed _ F248 = 2.0 kcal/ 
mol) relative to the open state (E open _ F2 48A — E open _ F248 = 
— 0.01 kcal/mol), an effect caused in part by the ab- 
sence in this case of robust stabilizing electrostatic inter- 
actions with the S5 and S6 transmembrane helices in 
the closed configuration. A dwell-time analysis showed 
in this regard that the mutation F248A caused a left dis- 
placement of the closed-time distribution curve with 
the main distribution centered at 30 ms for wild type 
shifted to 6 ms for F248A (Fig. S3) . These changes were 
coupled to a right shift of the open-time distribution 
from 7 ms for wild type to 50 ms for the F248A mutant 
(Fig. S3). However, the effects observed on the open- 
time distribution remains substantially less important 
than the changes seen with the F248W mutant. This 
conclusion can be extended to the F248V mutant. Col- 
lectively, the correlation illustrated in Fig. 7 (B and C) 
between ln(Keq_mutant/Keq_wild_type) and AAE oc 
tends to support a model whereby the pore helix is a 
major determinant in setting Pomax through non- 
bonded interactions involving F248 in the closed state. 
Of interest, a gating mechanism that takes into account 
possible interactions between the channel pore helix 
and the S5 transmembrane segment is in line with the 
observations reported for the hERGl potassium chan- 
nel where the interaction (hydrogen bonding) between 
H562 of S5, equivalent to W216 in KCa3.1, and the 
amino acids T218 (1244 pore helix in KCa3.1) and S221 



(T247 pore helix in KCa3.1) of the pore helix is func- 
tionally important for deactivation and Vi/ 2 of activation 
(Lees-Miller et al., 2009). It was also proposed that hy- 
drophobic interaction between L273 (L215 residue in 
S5-KCa 3.1) and V310 (F248 in KCa3.1) in the KCNQ1 
potassium channel stabilizes the open state of the selec- 
tivity filter, so that disruption of these interactions by 
mutating V310 to the small residues Ala or Gly could 
destabilize the pore open state, causing enhanced rate 
and extent of inactivation (Seebohm et al., 2005). How- 
ever, such behavior contrasts with the current proposal 
of F248 interacting with the S5 transmembrane helix 
leading to a stabilization of the channel closed state. 

Additional evidence in favor of wild-type Pomax being 
governed by F248 interacting with S5 in the closed state 
came from measuring the Pomax sensitivity to L211 mu- 
tations, as interaction energy calculations and the struc- 
tural model presented in Fig. 3 predicted a strong 
coupling energy between F248 and L211 for the 2A79- 
based KCa3.1 structure. Although a positive correlation 
was observed between the side-chain volume of the resi- 
due at position 211 in S5 and Pomax with a Pearson 
correlation coefficient of 0.75 ±0.15, our mutation study 
failed to show a strong perturbing effect of L211 muta- 
tion on Pomax. However, a significant increase in Pomax 
relative to wild type was observed with L21 IF, suggesting 
that L211F can come in proximity of F248 and stabi- 
lizes the channel open configuration through aromatic- 
aromatic interactions. Clearly, the wild-type KCa3. 1 Pomax 
is governed more by the interactions between F248 and 
W216 for the channel closed configuration than by the 
interactions between F248 and L211 for the channel in 
the open configuration. These interactions constitute 
potential targets for a pharmacological control of KCa3.1 
in accordance with the observation that mutating in 
KCa2.x channels the residue equivalent to L215 in S5 
affects the action of the KCa2.x activators GW542573X 
and CM-TPMF (Hougaard et al., 2012). However, an 
F248-based model cannot totally account for the lower 
Pomax reported for KCa3.1 compared with members of 
the KCa2.x channel family, as both KCa2.x channels 
and KCa3.1 have a Phe and a Trp at positions equivalent 
to 248 and 216 in KCa3. A lower Pomax could still arise 
from differences in the overall interaction between the 
pore helix and the S5 segment if one considers that 
the amino acid composition of the pore helix and S5 
transmembrane segment for KCa2.x and KCa3.1 differs 
by 46 and 70%, respectively, with four of the seven pore 
helix residues at the interface with S5 not conserved be- 
tween the KCa2.x and KCa3.1 channels. In addition, we 
showed in a previous work that Pomax could be in- 
creased by mutating residue S367 in the KCa3.1 CaM- 
binding domain to the equivalent Thr residue in KCa2.x 
channels (Morales et al., 2013). The difference in Pomax 
seen between KCa3.1 and KCa2.x channels likely reflects 
a global difference in the channel open configuration, 
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despite an identical CaM-based mechanism at the ori- 
gin of channel activation. 

Coupling F248 to the S6 transmembrane helix 
Although the proposed analysis highlights the impor- 
tance of F248 interacting with the S5 transmembrane 
helix in the closed state to Pomax, MD-based calcula- 
tions of nonbonded energy between F248 and residues 
in S6 showed strong attractive interactions with the Gly 
hinge residue at 274 and to a lesser extent with C277, 
M273, and T270 on S6 in both model structures. The 
strong interactions with G274 essentially reflect the 
tight packing between the F248 side chain and the G274 
backbone, suggesting a possible allosteric coupling be- 
tween the conformational changes induced by Ca 2+ at 
the level of the CaM-KCa3.1 complex and the gating 
events at the selectivity filter. Significant coupling en- 
ergy values were also obtained in double mutant cycle 
experiments using the F248W-M273W and F248W- 
C277W systems. The extent of the coupling depended 
on the size of the residue at 273 and 277, as no signifi- 
cant AAG values were obtained with the F248W-M273A 
and F248W-C277A double mutants. In fact, Pomax val- 
ues for the C277W and C277A (not depicted) mutants 
did not appear statistically different from wild type, thus 
ruling out perturbing interactions of the C277W or 
C277A (see Fig. 10) with F248. We noted, however, that 
the mutation C277W could reverse the Pomax increase 
resulting from the substitution of F248W, suggesting 
that aromatic-aromatic interactions could in this case 
stabilize the pore helix closed configuration (Fig. 11). 
MD simulations demonstrated in this regard that the 
mutation C277W on the KCa3.1-F248W template caused 
a greater stabilization of the closed (Eclosed_F248W_ 
C277W-Eclosed_F248W = -2.3 kcal/mol) than open 
(Eopen_F248W_C277W-Eopen_F248W= -1.3 kcal/mol) 
configuration for a AAEoc relative to F248W of +1.0 kcal/ 
mol. These simulations are thus in agreement with the 



observed decrease in Pomax coming from the mutation 
C277W on the KCa3.1-F248W template, an effect caused 
for the most part by an increase in electrostatic interac- 
tions between F248W and C277W (Fig. 11). However, 
the conclusions of our MD simulations need to be inter- 
preted taking into account that our models do not in- 
clude the transmembrane segments S1-S4 and thus 
provide a partial description only of the residues in the 
vicinity of C277. Although our observations confirmed 
that C277W is within the interaction range of F248W, 
we cannot conclude that there is an important contri- 
bution of C277 interacting with F248 in setting Pomax 
for KCa3.1 wild type. Our results also show an impor- 
tant increase in channel activity caused by the M273W 
substitution with a Pomax of 0.95 + 0.005 (n = 3) . These 
observations are in line with results obtained for the 
Kir3.4 (S143T) channel, where mutating to either a Val 
or Iso the residue LI 74 equivalent to M273 in KCa3.1 
resulted in constitutively active channels (Rosenhouse- 
Dantsker and Logothetis, 2006). On the basis of these 
observations, it was argued that the position above the 
central Gly could play the role of the gating hinge and 
control gating. 

Toward a global gating mechanisms for KCa3.1 
Our results added new elements to the complexity of 
the KCa3.1 gating mechanism by highlighting the im- 
portant role played by the interactions between the S5 
transmembrane segment and the channel pore helix in 
setting Pomax. In contrast to voltage-gated channels 
where gating originates from the movement of the S4 
transmembrane helix, KCa3.1 gating is initiated by the 
formation of a CaM-KCa3.1 complex in the C terminus 
directiy connected to the S6 transmembrane helix. The 
formation of the CaM-KCa3.1 complex most likely in- 
volves the formation of the dimeric structure between 
two adjacent KCa3. 1 subunits and potentially the forma- 
tion of a dimer of dimers, as discussed in a previous 
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Figure 1 1 . Structure representations of single 
and double mutants supporting close proxim- 
ity between C277 and M273 of S6 and F248 
of the pore helix. Representation of the chan- 
nel in the closed state (3BEH template). 
(A) The relative orientation of the Thr (W) for 
the F248W-C277W mutant suggests possible -ir-ir 
interactions that would stabilize the pore helix in 
the closed configuration, thus accounting for the 
mutation C277W preventing the Pomax increase 
coming from the F248W mutation. (B) Also il- 
lustrated is the tight packing coming from the 
mutation M273W, which confirms the proximity 
of F248 and M273, as demonstrated through the 
double mutant cycle analysis in Fig. 10. Illustra- 
tions by Discovery Studio Visualizer (Accelrys) . 
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work (Morales et al., 2013). Several studies on Kv chan- 
nels have demonstrated that interaction between the S6 
transmembrane segments and the S4-S5 linker at the 
cytoplasmic side of the membrane is essential for elec- 
tromechanical coupling (Lu et al., 2002; Ferrer et al., 
2006; Prole and Yellen, 2006; Choveau et al, 2011; 
Labro et al., 201 1) . This interaction would be crucial to 
couple the voltage-sensor movement in S4 with the acti- 
vation gate at the C-terminal end of the S6 transmem- 
brane segment (Long et al., 2005b; Labro et al., 2011). 
Given that the Ca' + -dependent CaM-KCa3.1 complex is 
connected to the S6 and not S4 transmembrane seg- 
ment, one would not expect a significant contribution 
of the S4-S5 linker region to the KCa3.1 gating process. 
However, the closed- and open-state models of KCa3.1 
illustrated in Fig. 1 (A and B) indicate that the S4— S5 
linker encircles the channel pore at the level of the bun- 
dle-crossing region, so that outward movements of the 
S6 transmembrane segment can be transmitted to the 
S4-S5 linker. A displacement of the S4— S5 linker could 
in turn produce a rotation and translation of the S5 
transmembrane helix and affect the interactions be- 
tween S5 and the channel pore helix. Our results clearly 
show that these interactions are crucial to the energy 
balance between the channel open and closed configu- 
rations. A possible gating process for KCa3. 1 would thus 
consist of the formation of a CaM-KCa3.1 complex in 
the C terminus leading to a displacement of the S6 
transmembrane helices that would induce a movement 
of the S4-S5 linkers, resulting in a rotation/displace- 
ment of the S5 transmembrane segments and to a modi- 
fication of the interactions between S5 and the pore 
helix. This process would work in parallel with a direct 
coupling of S6 with the pore helix and the selectivity 
filter region as documented in this work. This mecha- 
nism could explain how the conformational change in- 
duced by Ca 2+ binding to CaM can be transmitted to the 
selectivity filter region, insuring gating. Validation of 
this model would, however, require an in-depth evaluation 
of the role played by the S4— S5 linker in KCa3.1 gating. 

Conclusion 

This work demonstrates for the first time that aromatic- 
aromatic interactions between the pore helix and the 
S5 transmembrane segment in KCa3.1 control the chan- 
nel Pomax. Thus, interfering with these interactions 
represents a promising approach to modulate KCa3.1 
channel activity. Therefore, the interface formed by the 
pore helix and the S5 transmembrane segment repre- 
sents a valuable site that can be used for the design of 
KCa3.1 potentiators. 
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